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Abstract 

A method for treating ferromagnetic chains coupled with antiferromagnetic 
interactions on an hexagonal lattice is presented in this paper. The solution of 
the ID part of the problem is obtained by classical transfer-matrix while the 
coupling between the chains is processed by mean-field theory. This method 
is applied with success to the phase diagram and angular dependence of the 
critical field of CsNiFs. Results concerning the general influence of single- ion 
anisotropy on the magnetic ordering of such systems are also presented. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Mean-field theory has long been a useful approximation for the study of phase transi- 
tions in a wide variety of systems, especially magnetic ones.S This comes from the great 
simplification of dealing with an averaged system instead of explicitely taking into account 
each individual interaction. As is generally the case, the interest here in using mean-field 
theory is its ability to predict phase transitions and to follow the evolution of relevant quan- 
tities close to the transition point. Unfortunately, transition points and critical exponents 
extracted from the mean-field theory are found to be significantly different from the ex- 
perimental ones, especialy when the effective dimensionality is small. This is due to the 
averaging process involved which neglects the important critical fluctuations close to the 
transition point. 

Although mean-field theory has some problems when applied directly to low dimen- 
sionality systems, it can be a very good approximation if used in conjunction with other 
techniques. By definition, quasi-one-dimensional (quasi-lD) systems have a direction in 
which the energy scale is much larger than in the others. Some numerical techniques such 
as classical transfer-matrix,ii'iil quantum transfer-matrix,i and Bethe- AnsatJ can be used 
to solve almost exactly the ID part of the problem and so, instead of using mean-field for 
all the interactions, it can be restricted to the coupling between the chains.! By doing this, 
important fluctuations, although not the critical ones, will be included and a more accurate 
solution can be expected. 

The present work is mainly concerned with CsNiFa and, to some extent other equivalent 
systems. In this hexagonal insulator, the 5 = 1 Ni^"'' ions are ferromagnetically coupled 
in chains, along with the F~ ions. The chains are well separated from each other by large 
Cs+ ions and are coupled by antiferromagnetic interactions. This arrangment causes a 
large spacial anisotropy between the Ni ions, the ratio of the distance between these ions 
along the chains and in the basal plane being 0.4205, with an interchain separation of 
6.27 A. Anisotropy is thus present in the magnetic properties, the intrachain super-exchange 
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(J|l = 20 K)i being much larger than the interchain one. The exact ratio between the 
super-exchange interactions is difficult to estimate because of the strong dipolar interaction 
arising between the ferromagnetic chains. It is important to note that dipolar field is almost 
nonexistant if the ordering along the chains is antiferromagnetic (AF), as it is in many other 
ABX3 compouds such as CsMnBrs and CsNiCls.S 

This strong dipolar field is responsible for the particular planar AF arrangementS of 
CsNiFs which is different from the expected 120° structure of a system with only AF inter- 
actions on a triangular lattice. This magnetic phase occurs for temperatures smaller than 
T/v — 2.77 K.lli. Lussier et al^ have also shown a strong dependence of the critical field as 
a function of the angle between the magnetic field and the chain axis. As an example, at 
T = 2 K, Hc± — 0.27 T while for the other direction Hc\\ — 2.3 T. An unsuccesful attempt 
was made to fit this peculiar angular dependence of the critical field using a mean-field model 
similar to the one of Refs. |ll] and [1^. In this paper, it will be shown that the experimental 
angular dependence is typical of a system formed by XY ferromagnetic chains with AF 
coupling between them on a triangular lattice. 

The Hamiltonian describing the magnetic properties of such systems is given by: 

n = - j|| ■ s.+i -j± J2 Sa ■ s, - gi2Bii ■ E + ^ E(Si)' + (i) 

i <a,b> i i 

The ffist, third and fourth sums in equation are along the chain while the second one 
is between nearest neighbors in the plane perpendicular to the chains axis. The strong XY 
behavior of CsNiFs comes from its large single-ion anisotropy (D) of 8.5 K.i The g-value 
was set to 2.2. 

As mentioned previously, the usual mean-field theory uses averaged quantities for the 



three lattice dimensions. An example of this is found in Ref. 0. The goal of the present 
work is to develop a method which treats the chains by classical tranfer-matrix and uses 
mean-field for the coupling between them. This method, refered to as MFTM (Mean-Field 
Transfer-Matrix), is presented in section II. It is followed, in Section III, by a brief review 
of the transfer-matrix algorithm used to calculate the magnetization and the susceptibility 



of a ID chain as a function of temperature and magnetic field. In Section IV, MFTM is 
applied to CsNiFa for the calculation of its magnetic-field - temperature phase diagram and 
angular dependence of the critical field. Some other general results, concerning single-ion 
anisotropy, are presented in Section V. 

II. MEAN-FIELD COUPLING OF QUASI-ID CHAINS 

Before presenting the calculated phase diagram of CsNiFs, it is useful to understand 
qualitatively the underlying mecanism of the magnetic order in this system. For simplicity, 
consider a starting point below Tc with a magnetic field oriented in the XY plane, large 
enough to destroy the AF order. In this phase, the paramagnetic one, the free-energy is 
dominated by a Zeeman term, — m ■ H, and the magnetization m is parallel to the applied 
field. As the field amplitude is lowered, this term is diminished up to a point where a 
perpendicular exchange term, proportional to —J±m'^, becomes equal to it. At this field, 
since < 0, an AF order develops and a finite angle between the magnetic sub-lattices 
appears. The critical field. He, is defined where this angle is equal to zero. It is important 
to notice the particular planar structure of the magnetic order in CsNiFs (Fig. [l|) which is 
different from the 120° structure of CsMnBra, an XY system having strong AF coupling 
along the chains. This difference is due to the relatively strong dipolar field originating 
from the neighboring ferromagnetic chains0. This dipolar field depends strongly on the 
angle between the magnetic sub-lattices and it can be shown quite easily that there is no 
interchain dipolar field present on a given site at He, where the magnetization of all the 
sub-lattices are parallel to the field. So, for small angles between the sub-lattices, the order 
starts to develop like a 120° system (see Fig. 0) but at a given angle, the dipolar field coming 
from the other chains becomes large enough to flip one of the sub-lattices according to the 
planar arrangement of CsNiFa (Fig. 0). 

The starting point of the theoretical description of such a system is the full three- 
dimensional {3D) Hamiltonian of equation (|I]). Since this Hamiltonian, applied to a tri- 
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angular lattice, involves up to three magnetic sub-lattices, it is a highly non-trivial task to 
solve it even numerically. On the contrary, the one-dimensional {ID) part of the Hamilto- 
nian consisting of first, third and fourth term can be calculated numerically quite easily and 
with good precision by the tranfer-matrix technique. The utility of mean-field theory for 
this problem comes from the fact that the ID part has a much larger contribution to the 
free-energy than the 3D one (| J||| ^ | Jx|).i Thermal averages of the ID Hamiltonian can 
be used to approximate the free-energy of the 3D one. Formally, this type of free-energy is 
obtained using the Bogoliubov inequality:0 

^ < -^T = -^ID + - nw),D ■ (2) 

With the three magnetic sub-lattices of Fig. |^, and with a proper counting of each link 
in the J± term, this gives the following trial free-energy: 

i 

(3) 

i 

The thermal averages J2i{^ia) are just the magnetisation m of the sublattice a. Using the 
angle definition of Fig. ^ for 6 and defining as the angle between the field and the XY 
plane, the trial free-energy can be written as: 

J-'t = J-'iD — 3 J_L (2 cos 9 + cos 26)w? 

- {1 + 2 cos 6) cos (j)Hm± , (4) 

where ■m± is the component of m in the XY plane. 

It is also important to note that equation (^) is valid only for small values of the angle 6 
where the dipolar field has no infiuence on the order. Minimising (^) with respect to 6 gives 
sin 6' = or, 



cos^^ 



1 vcos (f)Hm^ 



rn? 



+ J 



(5) 



2J_ 

This results is reasonable. At high fields, ^ = 0, indicating the field-induced ferromag- 
netic (paramagnetic) phase. In zero field, and for a sufficiently low temperature, one gets 



the 120° configuration. Of course, tliis type of structure is not applicable for CsNiFs but 
it reflects the ordering mecanism which is valid only for small values of 6. As indicated 
above, He is defined where cos 6 = 1 and it can be obtained by solving self-consistently the 
following equation: 

— 3J^m^ = cos (j)Hm^ (6) 

Each side of equation (j^) represents a contribution to the free-energy JF (to within a 
minus sign). The right-hand side is the magnetic field contribution — m ■ H, the dominant 
one in the paramagnetic phase, while the left-hand is the 3D contribution of the ordered 
phase. Since m is bounded to unity, the field contribution will eventually dominate at high 
field and the phase with the lowest free-energy will be the paramagnetic one. The strength 
of the ordered phase is controlled by the parameter J±. In some systems, like CsNiFs, 
this parameter is not well known but it can be evaluated with this technique by fitting to 
experimental data. 

From the limiting behavior of equation (^) at low field, it is also possible to calculate, 
within MFTM, the transition temperature in zero fieldjl T^. This gives: 

XlD(Te,0) = ^, (7) 

where Xid is the ID susceptibility. An interesting application of equation (|^) is the possi- 
bility to calculate directly the value of J± needed for a given T^. In the presents work, this 
is the only procedure used to set J±. 

The last item that needs to be solved in this mean-field theory is the value and the ori- 
entation of the external magnetic field Hext- Up to now, the field included in the calculation 
was the effective field applied to ID chain. This is not the true external field; one needs to 
add a dipolar contribution Hdip and a 3D contribution, H^d, coming from the influence of 
the neighboring chains. Using Hid to designate the field used previously in this paper this 
gives: 

Hext = Hi£) — H3/) — Hdip ■ (8) 
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Equation (§) is shown graphically in Fig. |^. This figure also shows two new angle definitions 
between fields and the XY plane: O^xt for Hext and 6^ for H^d and iidip- The evaluation 
of is straighforward it is just the number of near neighbors times J±. It gives: 

H3D = 6Jxm (9) 

The dipolar field Hdip can be evaluated easily, if the magnetization of all the sub-lattices are 
parallel, since there is no interchain contribution. For m in the plane perpendicular to the 
chain this gives: 

Hdip = -32.3/^m , (10) 

with 

(11) 

Since CsNiFs has a quite strong XY caracter, the component of m in the XY plane was 
always by far the greatest for all the angles used in the calculation so that the numerical 
constant used in equation (p!oD was good to a few percent in all cases. 

The fields Hsij and 'Rdip are important in evaluating the external field tlext but they 
are not included in the calculation of by equation (^). The first one, Hs^i, is implicitely 
accounted for in the left-hand side of equation (|^) while the second one, H^jp, obviously 
turns with m so there is no free-energy changes due to it. 

The phase transitions are obtained by the solution of equation (^ which implies a prior 
knowledge of m.{T,H), the ID magnetization. An effective and flexible way of calculating 
m is the transfer-matrix technique, which depends on the various ID parameters like the 
single-ion anisotropy (D) and the parallel exchange interaction {J\\)- 

III. TRANSFER-MATRIX TECHNIQUE 

Since the present work relies heavily on the tranfer-matrix technique, a brief summary 
of this method is presented here. This study follows the work of Blume et a/.i and assumes 
that the classical spin Hamiltonian for an spins chain of magnitude 5* can be written as: 



with 



and 



N 



i=l 



J 



J\\S{S + 1) rj^l^p 

\J\\S{S + 1)\ lj=-l^AF 



(12) 



(13) 



(14) 



D 



DS{S+1) rD <0^ Ising 



In 



|J||5(5+1)| ^D>{)^XY 



(15) 



H 



IJII^(^ + 1)I 



(16) 



/3 = /5|J||5(5 + l)| 



(17) 



The unit vectors 5*, are oriented along the spin direction. Using these definitions, the 
partition function can be written^ as a product of Boltzmann factors, 



N 



Z = ^exp[/35:\/(5i,5, 



i+lj 



i=l 



N 



= Enexp[/3m,5m) 

{si} «=1 

Since V is translationnaly invariant, the term in the product does not depend on the site 
chosen and Z is given by: 



Z = TtK 



N 



(18) 



with 



(19) 



Since is a very large number, only the largest eigenvalue Aq of K will contribute signifi- 
cantly to Z. 
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For classical spins, the transfer-operator K needs to be mapped onto some discretisation 
of the spherical coordinates in order to evaluate it numerically. Following the procedure 
found in Ref. ^, which is suitable for broken rotational symmetry, one finds the ID magne- 
tization to given by: 

m = E E W!>Wf^o{e,, (j),)^o{0^, 0,) [cos{e,)z + sm{e^) cos(0,)x + sin(^,) sin(0,)y], (20) 

i=l j=i 

where \l/o is the eigenvector of K corresponding to the largest eigenvalue Aq and, Ne, N^, 
W^, Wf, are respectively the number of coordinates and the Legendre weights of the discret 
spherical coordinates 6i and 0j. In the absence of a magnetic field, or when the field is 
along the chain axis, a simplification using the rotational invariance of V can be usedE 
The reduced number of coordinates for K greatly improves the amount of computing time 
needed for a given precision. Such a simplification can be used to evaluate the susceptibility 
X of the zero field equation (0). Using the results found in Ref. |]the ID susceptibility can 
written as: 

xiL (g) = ^ E ,2 '° ^,2 [ E cos(^.) vi/oo {e.)%o{e^)] 

Aqo — ^AooA/o cosg + A;o j=i 
XtU^l) = f E ,2 2A A cL^X^ E^^cos(g0^oo(^O^u(g.) . (21) 
When using the rotational invariance as a simplification, one needs to introduce a second 
index in order to specify the azimuthal dependence. This second index is implicitely taken 
into account with the broken rotational symmetry algorithm. Typical matrix sizes are 
28^ X 28^ for the broken rotational symmetry algorithm while they are only 64 x 64 for the 
rotational invariance one. In both cases the amount of computing is not excessive, ranging 
between seconds and minutes per (T, H) point. Of course this figure increases rapidly with 
matrix size. 



IV. APPLICATION TO CsNiFg 



The main goal of this work is to calculate the magnetic phase diagram of CsNiFs . Figure ^ 
presents a comparison between the experimental phase diagram obtained with ultrasound 
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by Lussier et a/.0 and the phase diagram calculated with the MFTM method. In this 
figure, the circles are the original experimental data and the squares are from the MFTM. 
The diamond at zero field is set to 2.77 K by an appropriate choice of J± according to 
equation (|^. The value found for J± using this procedure is 0.0253 K. For this phase 
diagram and all the subsequent MFTM data on CsNiFs it will be the only value of J± used. 
The overall agreement between theory and experiment is good considering the fact that this 
mean-field theory neglects the 3D fiuctuations which are not included in the transfer-matrix 
calculation and are the relevant ones close to Tc. It is quite obvious from the development of 
MFTM that He scaled as m and so He ~ (Tn — T)^. The absence of the 3D fiuctuations thus 
refiects itself in the critical exponent (3, giving the usual mean-field value of 1/2 value for 
the theoretical phase diagram. This compares with an experimental estimateB of (3 = 0.31. 

The large XY single-ion anisotropy in CsNiFa has an important effect on the angular 
dependence of He. Crudely, one can say that when increasing the field for T <^ the 
transition occurs when m saturates so that the left side of equation (^) becomes bounded 
and the H term at the right can dominate. Because of the XY anisotropy of CsNiFs, 
m will grow rapidly along the projection of H in the XY plane and it will saturate for 
approximately the same field for small angles of the field out of the plane. For large angles, 
the projection of H in the XY plane will be small and m will grow more slowly giving a 
larger He- This type angular dependence is shown in Fig. |^. In this figure, the experimental 



data of Ref. |T2| and the MFTM results are represented respectively by the circles and the 
diamond. The agreement, already quite good between these data, becomes excellent when a 
possible correction for a misalignement of the experimental plane of rotation of 6° is included 
(circles). This correction, although large, is not unreasonable. One should also consider the 
nature of the theoretical method used as a source of error, especially close to = 90, a 
pathological angle for the present method since the coupling between and H then tends 
to zero. As mentioned earlier, even for the largest angle of the magnetic field out of the XY 
plane used, m was always out the plane by less than 10° so that the numerical constant in 
equation (p!oD was valid to a few percent. The good agreement found in Fig. || constrasts 
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sharply with the failure of full 3D mean-field theory, as mentioned in Ref. 



12 



V. SINGLE-ION ANISOTROPY 

More generally, in the absence of a magnetic field, the paramagnetic phase transitions 
are determined by the magnetic susceptibility of the system according to equation (^. 
When a single-ion anisotropy is present, like in the Ising and XY cases, the susceptibility 
is not isotropic and only the largest component needs to be taken into account. For an 
Ising system it is the longitudinal component Xihi^) while for an XY one, it is one of the 
transverse components Xihi^) The equation (|^) is particularily useful since from 

the knowledge of the susceptibility at a given temperature one can extract the perpendicular 
exchange interaction J± needed for a phase transition at this temperature. Such plots, for 
various values of D in a system having J|| = 1 and 5* = 1 are shown in Fig. ^. One 
can observe in this figure three limiting cases; the Heisenberg one at high values of J±, 
the XY one at intermediate values and the Ising one at low values. This can be easily 
understood, the general effect of the single-ion anisotropy is to favor the order by increasing 
the susceptibility when reducing the degrees of freedom of the spin system. It is especially 
true if the anisotropy is of the Ising type since the degrees of freedom become discrete. 
The XY anistropy has a similar influence but it is less pronounced, the degrees of freedom 
being still continuous. Of course, at high enough temperature (T ~ D/ks) the Heisenberg 
behavior is recovered in all cases. 

There is a functional difference between the continuous and discrete degrees of freedom. 
For the continuous ones, the low J_|_ behavior is characterised by a power law of the from 
Tc ~ J^^. The dotted hne, having a 1/2 slope, helps appreciate this. On the contrary, the 
discrete one (Ising) shows a more complex exponential form. These results are in agreement 
with the arguments of Villain et a/.0 Their result for the Ising case is shown by the dashed 
curve. 

The last aspect investigated in this paper is the influence of the amplitude of the XY 
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single-ion anisotropy on the whole phase diagram. As seen in the previous paragraph an 
increase of the anisotropy will enhance the magnetic order as a function of the temperature 
by increasing the susceptibility. As a function of the magnetic field, for T ^ Tc, the 
anisotropy should have a much smaller effect since it is not implied directly in the phase 
transition. Figure |^ shows such phase diagrams for a magnetic field in the XY plane and 
for three values of D; 0.1, 0.2 and 0.5. The value of J± has been set arbitrarily and the 
same lattice parameters than CsNiFs have been used. In zero field, the predicted behavior 
is observed, the highest Tc being with D = 0.5 while the lowest is with D = 0.1. For a 
constant temperature smaller than Tc (like T /2J\\ = 0.06), one can also observe the small 
effect of D on He- 



VI. CONCLUSION 

In this work, a method for treating antiferromagnetically coupled ferromagnetic chains 
by mean-field theory has been developed for a system with a hexagonal lattice. The solution 
of the ID part of the problem has been obtained by a generalised tranfer-matrix algorithm 
which is suitable for broken rotational symmetry.! Although special care has been needed 
to correctly account for the large dipolar field originating from the ferromagnetic chains, it 
has been possible to reproduce with succes the phase diagram of CsNiFs by only adjusting 
the value of J±. The peculiar dependence of the critical field of CsNiFa as a function of the 
angle between the magnetic field and the chain axis has also been reproduced successfully,^ 
in sharp contrast with 3D mean-field theory. 

More generally, the dependence of the critical temperature as a function of J± has been 
established for various values of the single-ion anisotropy D. For the cases having continuous 

1/2 

degrees of freedom, like the XY and Heisenberg ones, a Tc ~ J/ has been found while for 
the Ising case a more complex form has been obtained. These results are all in agreement 
with the ones obtained by Villain et al in Ref . |1^. For the XY case, the effect of the size of the 
single-ion anisotropy on the overall shape of the phase diagram have also been investiguated. 
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The results are in agreement with some simple physical arguments, an increase of Tc with 
D and a rather small effect on H^. for temperatures much lower than T^. 

On the basis of the success of this method for the phase diagram of CsNiFa, some exten- 
sions can be considered. The phase diagram of AF-chains compounds with XY behavior, 
like CsMnBr3,0 could be examined although some complications with the transfer- matrix 
algorithm, coming from the difference in the effective field of the two staggered sub-lattices, 
are expected. Since successive phase transitions occur, the present method, based on linear 
response theory, must be generalized so that the free-energy itself is calculated.0 A more 
challenging prospect is the application of this method to easy-axis AF chains systems like 
CsNiCl,.§ In this case, the single-ion anisotropy plays an active role in the magnetic order- 
ing so it cannot be taken into account only by the transfer-matrix solution of the ID part 
of the problem; it must also be included in the calculation 3D free-energy. This will add at 
least an order of magnitude to the complexity of the problem. 
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FIG. 1. Planar arrangement of the two magnetic sub-lattices of CsNiF3 for = and T < Tc 
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FIG. 2. Ordering of a three antiferromagnetic sub-lattice system at a magnetic field close to 
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FIG. 4. Comparison between the experimental phase diagram of CsNiFs given in Ref. ^ (cir- 
cles) and the one calculated by the method described in this paper (squares). The diamond 
represent the zero field point used to set J± according to eqn. 



17 



♦ Theory 

O Experiment 

□ Experiment with correction 



□ o 
□ o 



15 30 45 60 75 90 
Angle (degrees) 



FIG. 5. Comparison between the experimental angular dependence of He at T = 2 K given 
in Ref. |l^ (circles and squares) and the one calculated by the method described in this paper 
(diamond). The squares represent the original data (circles) corrected for a plane of rotation 
misalignment of 6°. 
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FIG. 7. Magnetic phase diagrams for various D values; 0.5 (circles), 0.2 (squares) and 0.1 
(diamond). The J± has been set arbitrarily. 
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